4 i heugceed wee. 


TID-4500, UC-53 
Meteorology 


I 


LAWRENCE LIVERMORE LABORATORY 
University of California Livermore, Caliomia/94550 


UCRL--51328 


AN INVESTIGATION OF SCAVENGING OF RADIOACTIVITY 
FROM NUCLEAR DEBRIS CLOUDS: RESEARCH IN PROGRESS 


W. K. Crandall 
C. R. Molenkamp 
A. L. Williams 
M. M, Fulk 
R, Lange 
J. B. Knox 


MS, date: January 17, 1973 


NOTICE 
This report was prepared as an account of work 
sponsored by the United States Government, Neither 
the United States nor the United States Atomic Enetgy 
Commission, nor any of their employees, nor any of 
theit contractors, subcontractors, or their employees, 
makes any warranty, express of implied, or assumes any 
legal liability or responsibility for the accuracy, com- 


Pleteness cr usefulness of uny information, apparatus, 
product of process disclosed, of represents that its use 


woutd not infringe privately owned rights, 


pisreemireiny ne TUTE FeONTARES SS se 


¢ . 
l tt 


a 


Ww 


Contents 


Abstract 
Introduction 


Summary of the Pertinent Findings from Potential Exposures frum 
Low-Yield Free Air Bursts F ; 5 ; ; ; ; 


Dynamic Effects of Convection a . . : : . . 
A Microphysical Description of Piecipliation Scavenging . . 
Thoughts on Altering the Scavenging of Debris Particles : 

Dry Deposition . ‘ F . ; 7 ' . 


A Three-Dimensional seas ey ae Diffusion Particle-In-Cell 
Code (ADPIC) . 3 A i 


Future Work 
Acknowledgments , 
References : ‘ ‘ , 


-iii- 


Saha 


AN INVESTIGATION OF SCAVENGING OF 
RADIOACTIVITY FROM NUCLEAR DEBRIS CLOUDS: 
RESEARCH IN PROGRESS 
Abstract 


The publication Potential Exposures 
from Low- Yield Free Air Bursts, 
addressed the problem of rapid deposition 
of radioactivity from nucleai debris 
clouds at various distances from ground 


1) A study of the dynamics of cumulus 
cloud interaction with a nuclear 
debris cioud to determine possible 
enhancement or reduction in peak 


surface concentrations. 


zero, The results showed that large 2) A study of the microphysical inter- 
values of infinite exposure and large actions of debris particles with 
exposure rates could occur at the surface cloud droplets and raindrops to 
under the following special circumstances: determine rainout and washcut 
the radioactivity is distributed in the efficiencies, 
lower troposphere, it is scavenged at an 3) A study to determine the feasibility 
early time by a rain-bearing system and of rendering the debris particles 
rapidly deposited over the ground, and it non-wettable, thereby effectively 
is confined to a relatively small area, making the rainout coefficient 
More recent studies have been and are negligibly small, 
being made at LLL to make it possible 4) A study of the potential effects of 
to assess 1) the reasonability of the re- dry deposition due to the detonation 
sults of this earlier publication by con- of larger numbers of low yield 
sidering in greater detail the physical devices within a confined area, 
processes involved and 2) the extent to 5) The development of a three- 
which the results may be changed by man- dimensional atmospheric diffusive 
made modifications. This paper is a prog- particle-in-cell code capable of zal- 
ress report on the research projects cur~ culating the time-dependent distri- 
rently underway at LLL which are directed bution of air pollutants under con- 
toward quantifiable assessinents for client ditions of turbulent diffusion and 
needs, Specifically, these projects are: wind shear, 

Introduction 


A major conclusion from the report 
Potential Exposure from Low-Yield Free 
Air Burst was that the hazard of close- 
in exposure from these burats may be 
appreciable when scavenging by precipita- 


ite 


tion occurs, For example, rainout of the 
radioactivity within the debris cloud of a 
1-kt nuclear detonation at 100 km dcewn- 
stream of ground zero may result in 


external gamma dose rates as high as 


100 R/nr and infinite whole body doses of 
1000 R orhigher, A study by Conner and 
Phillips” further substantiates the hazard 
inherent in these low-yield bursts. 

From the interest generated by these 
papers, it has become apparent that a 
greater effurt is required to obtain more 
detailed and more accurate assessments 
of the hazards involved. With this need 
in mind, the Atmospheric Sciences Group 
at Lawrence Livermore Laboratory some 
months ago initiated work on a suite of 
vesearch projects to iraprove understand- 
ing and assessment capabilities pertain- 
ing to the scavenging of debris from 
nuclear clouds, It is the intent of this 
paper to present an outline of eaeh of 
these projects, their current status and 
some preliminary results, 

The specific research problems to be 
discussed are as follows: 

(1} & model of the interaction between 
a curn.us cloud and a nuclear debris 
cloud to determine possible enhancement 
or reduction in surface concentrations 
due to relative positions and motions, 
Among the assumptions used in the early 
ealculations was that the vertical integral 
through cloud-center was deposited on the 
ground and integrated exposures and 
exposure rat..s were calculated based 
upon this deposition, It is possible, how- 
ever, that debris will be entrained into 
a developing cumulus system and that the 
resultant vertical integral of the debris 
concentration in the cumulus cloud may 
be enhanced by this redistribution of 
radioactivity. Conversely, the resultant 
vertical irtegral may be decreased de- 
pending upon the relative vertical positions 
of the scavenging system and the debris 


cloud, 


“Qe 


(2) A study of the microphysical inter- 
actions of debris particles with cloud 
droplets and raindrops to determine 
It is known that the 


rate at which the debris is removed from 


deposition rates. 


a cloud depends upon the rain rate, the 
spectrum of particle and raindrop sizes, 
the relative fall velocities of the raindrops 
and the particles, the efficiencies of 
capture for a given particle size, and the 
chemical properties of the particles. 
Frem these parameters washout and 


rainout efficiencies can be calculated 
and some measure of the effective- 


ness of the rainout process can be 
determined, 

(3) A study to determine the feasibility 
of rendering the debris particles non- 
wettable, 
air bursts appear to be insoluble and 
wettable, This being the case, particles 
with radii as small as 0.1 wm can serve 
Acccrding to 


Particles resulting from free 


as condensation nuclei, 
arguments by Fletcher® and experimental 
results by Twomey, particles whose 
surfaces depart from being completely 
wettable cannot serve as condensation 
nuclei under typical cloud supersaturations, 
Since nucleation scavenging appears to be 
the dominant rainout mechanism, a slight 
-eduction in t' e particle wettability could 
greatly reduce the potential hazard from 
rainout, 

(4) A study of the potential effects =f 
dry deposition (deposition caused by the 
dispersion of debris particles to the sur- 
face) due to the detonation of a large 
number of low yield devices within a con- 
fined area, 

(5) The development of a three- 
dimensional Atmospheric Diffusive 
Particle-in-Cell Code capable of 


calculating the time-depending distri- 
bution of air pullutants under conditions 
of turbulent diffusion and wind shear. 

A summary Of the pertinent findings of 


the report Potential Exposures from Low- 
Yield Free Air Bursts 
the benefit of readers not familiar with 


is included for 


this previous report, 


Summary of the Pertinent Findings from 


Potential Exposures from Low-Yield Free Air Bursts 


The risk due to fallout from a free air 
burst has generally been considered 
negligi’sle since few particles large 
enough to settle to the ground are created. 
Early time rainout or wasuout of radio-~ 
active debris from air bursts has been 
considered a potential hazard® but since 
choice of meteorological ccnditions has 
always been available during testing, the 
probability of damage from this phenomenon 
has been consiuered small or non-existent. 
If, however, it becomes necessary to 
detonate many low yield devices within a 
Smail area, at a time when choice of 
meteorology is not an option, then hazards 
due to rainout or washout become much 
more probable. 

Consider the following assumptions 
relating to the detonation of an all fission 
nuclear device: 

1. It is detonated at a low altitude but 
high encugh that there is no inter- 
action of the fireball with the sur- 
face (free air burst),° 
The resulting debris cloud rises to 
its stabilization altitude and is 
approximated in form by a right 
circular cylinder, 

The distribution of radioactivity in 
the debris cloud at stabilization is 
assumed Gaussian in both the radial 
and vertical directions, 

The debris cloud moves downstream 
at the speed of the wind at the 


w 
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stabilization altitude, 
The debris cloud undergoes no 


a 
. 


shear, 
The debris cloud is dispersed 


n 


radially and vertically by diffusive 
eddies in the atmosphere. 

7, At some point downstream, the 
debris in the debris cloud is sub- 
jected to scavenging by a rain 
bearing system, th-reby causing 
all or a portion of the radioactivity 
to be deposited upon the ground. 
This deposition process is assumed 
to be instantaneous, 

8, Once the radioactivity is deposited 

on the ground, it remains in place, 
Assuming the above sequence of events 

occurs, it is possible to calculate a rate of 
exposure and an infinite exposure to a per- 

son standing at the deposition point once a 

means of determining the areal concen- 

tration of the radioactivity is specified. 
For this study, the areal concentra- 
tion is defined as the radioactivity re- 
moved from a column of unit cross 
sectional area through the center of the 
debris cloud (also often referred to here 
as the vertical integral’). Thus as the 


“Depending upon use, the term vertical 
integral may refer to the radioactivity 
deposited, the dose rate due to this radio- 
activity, or the infinite dose due ta this 
radioactivity, and, unless otherwise 
specified, it refers to cloud center, 


fa jail Riese ada: Re 


debris cloud moves downstream, its 
vertical integral can be calculated as a 
function of time from the diffusion param- 
eters, and the pertinent exposures and 
rates can then be calculated, The vehicle 
for calculating the vertical integral is the 
Large Cloud Diffusion Code 2BPUFF.® 
Table 1 represents the input param-~ 
eters to 2BPUFF for vertical integral 
calculations of debris clouds from one, 
ten, and one-hundred kt detonations. The 
cloud dimensions, height or rise,’ and 
meteorology are based upon mean con- 
ditions. Two diffusion conditions, 
atmospheric dissipation (e€) and vertical 
diffusivity (K), are assumed in order to 
obtain a reasonable range cf values of 


the vertical integral as a function of 


integrals (converted to time of arrival 
dose rates) in R/fnr for the three yields, 
and Fig, 2 shows the infinite whole body 
exposures. The time-dependence has 
been converted to dependence on distance 
from ground zero in both figures, I[t is 
seen that extremely large exposures and 
rates can occur at distances even as far 
as 1000 to 2000 km for the 100 kt case, 
However, another importarit factor is that 
the debris clouds for the large~-yield 
bursts will rise to greater heights and 
some or all of the debris wiil no longer 
he available for scavenging because no 
rain bearing systems will reach high 
enough, Thus, based upon an estimate 
of types of scavenging systems, it has 
been assumed in this report that 100% of 


time.©9-19 Figure 1 shows the vertical the 1-kt debris, the lower 10% of the 
Table 1, Input cloud parameters for exposure calculations. 
Atmospheric 
dissipation, Vert cal diffusivity 
Center Thick- 2 Grose” 131 wind, —¢(ergs/e-sec) itor fuec) 
Yield height ness? Radius fiasicn t speed Siow Fast Stax Fast 
(kt) (m) — (m) (m) (pCi) {pCi} (km/hr) diffusion ‘diffusion diffusion diffusion 
1 «2,840 += -1760 920 4.44107 1.4x10!7 39.8 0.5 3 1000 10,000 
10 7,000 3060 . 2400 444x107! a.4xipl® 70.2 1 5 1000 10,000 
100 11,700 5340 000 4.44x 107 z4x10'? 72.9 0.7 


8corresponds to two standard deviations from sloud-center concentration. 


Patter H+} hr. 


3.5 1000 10,000 


Table 2, Potential surface exposure and downstream exposure rates. 


Infinite External gamma 
Fraction Distance whole-body exposure rate at 
Yield deposited downwind external gamma dose time of arrival 
(kt) (7%) (km) (rem) (R for) 
1 100 10 25,000 25,000 
100 200-1200 35-100 
1000 0.3~1.0 <0,01 
10 10 10 5000 7000 
100 1000-1500 150-200 
1000 1-4 0.015-0.05 
100 J 10 800 1200 
100 300-350 45-55 
1000 1,5-8 0.02-0.06 
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Initia] whole-body exposure rate — R/hr 


10 


10 


107! 
10 100 1000 
Distance from ground zero — km 
Fig. 1. Initial gamma exposure rate due to the deposition of the vertical integral at 
various distances from ground zero, The upper curve for each yield repre- 


sents the case of slow horizontal diffusion; the Lower curve represents the 
case of fast diffusion, 
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Infinite whole-body exposure — R 


Fig. 2. 


ee a 
100 


Distance from ground zero — km 


Vertical integral (infinite whole-body exposure) due to gross gammc radiation 
as a function of distance from ground zero. The upper curve for each yield 
represents the case of slow horizontal diffusion; the lower curve represents 
the case of fast diffusion, 
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10-kt debris, and the lower 1% of the 
100-kt debris are cvailable for scav- 
enging. 

With these limits as a basis, Table 2 
summarizes exposures and exposure 
rates for each of the three yields for 
distances of 10, 100, and 1000 km down-~ 
stream from ground zero, Observe that 
while the total exposures due to the 10- 
and 100-kt clouds are reduced by factors 
of 10 and 100 respectively, rather high 
values can still occur even out to 100 km, 


Dynamic Effects 


In the preceding section of this report 
it has been assumed that nearly all the 
debris in a vertical column from the 
ground up to the top uf a rain cloud is 
But the 


presence of rain clouds implies the 


scavenged hv precipitation. 


existence of dynamic effects which can 
greatly alter the vertical distribution of 
debris from that in the cloud-free environ- 
ment, 

Consider, for example, a simplified 
picture of the dynamics of a growing 
cumulus cloud, Moist unsaturated 
environmental air is drawn up through the 
base of the cloud and is accelerated up- 
ward because of buoyancy caused by the 
relvase of latent heat of vaporization, 
This indrawn air either remains a part 
of the young growing cloud or it leaves 
through the anvil top of a large fully 
developed cloud, The debris tends to 
move with the air so that the concentration 
of debris at all heights in a cumulus cloud 
developing in this manner would tend to 
be the same as the debris concentration 
of the environment in the vicinity of the 
base of the rain cloud, Now if the debris 


It should be recelled, however, that the 
10% and 1% availability of the 10- and 
100-kt cloud debris was based upon a 
specific case and may ke greater or less 
depending upon individual situations. For 
example, Conner and Phillips” ina 
Similar calculation founa that 100% of 
the radioactivity from the 10-kt cloud 
was available for scavenging based upon 
his calcutations of stabilization aeipht, 
cloud thickness, and height of scaveuging 


systems, 
of Convecticn 


concentration is maximum near the rain 
cloud base, the amount of debris in a 
vertical column in the rain cloud will be 
much larger than the amount of debris in 
a corresponding column of the clear 
atmosphere, When rain falis, moe than 
one times the vertical intergral of the 
debris concentration in the clear atmos- 
phere can be brought to the ground 
beneath the debris-cloud center, 1n addi- 
tion, the mass of air that passes through 
the rain cloud may be much larger than 
the mass of air in the cloud itself, further 
increasing the amount of debris that can 
be collected by drops and droplets and 
eventually brought to the ground, On the 
other hand, if the debris is primarily 
above the base of the rain cloud, it can 

be imagined that the cumulus cloud would 
grow up through the contaminated regions, 
pushing the debris aside and thereby 
preventing collection of debris by rain- 
drops, From these simple examples it, 
can be seen that the dynamics of individual 
convective cells are a maior factor in con- 
trolling the amount of debris that enters 


the vain cloud and interacts with drogs. 


As debris moves through a rain cloud, 
some of the radioactivity is collected by 
raindrops and cloud droplets, Some of 
the cloud droplets are collected by falling 
raindrops wile others evaporate when 
the cloud dissipates. The contaminated 
raindraps eventually fall to the ground, 
carrying with them ail of the radioactive 
debris taey have collected. Therefore, 
the debris in a rain cloud that is collected 
by rafiudrops is scavenged, and the radio- 
aciivity that does not get into raindrops is 
left in the atmosphere, 
microphysical processes by which debris 


particles get intc raindrops are the major 


Consequently the 


factors that determine what fraction of 
the debris passing through a precipitating 
system is actually scavenged. 

The remainder of this section is devoted 
to a description of a simple mcdel of 
cumulus convection that has been used to 
esti:nate the typical range of values for 
the ratio of the vertical integral of the 
debris concentration inside a rain cloud 
tc a Similar vertical integral in the cloud- 
free environment. This ratio of vertical 
integrals is given the name enhancement 
factor, although enhancement factors 
less than 1 correspond to cases in which 
the rain cloud contains less debris than 
the environment. In following sections 
of this report, the scavengability of debris 
varticiles and the microphysical processes 
by which debris is collected and brought 
to the ground with falling raindrops will 


be discussed. 


DESCRIPTION OF THE MODEL 


A one-dimensional steady-state cumulus 
dynamics model has been used to estimate 
the importance of the dynamics of cumulus 
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convection on the scavenging of debris 
from the atmosphere. This model was 
developed by Weinstein and Davis (1968).!! 
It follows a sample parce! of air rising 
through a rain cloud; the cloud as a whole 
is assuned to be made up of many of 
these parcels one on top of the other and 
each with the properties of the sample 
parcel at the appropriate level; hence the 
structure of the rain cloud does not change 
with time. The sample parce! is initially 
at cloud base and has a small upward 
velocity. The surrounding environment 
has a vertical profile of temperature and 
moisture content determined from a 
Calcula- 
tions of the cloud tempe ature, vertical 


standard radiosonde sounding. 


velocity, mixing ratio, cloud water content, 
hydrometeor water content, and undraft 
radius are made for the sample parcel at 
each level in the model cumulus cloud, 

The sample parcel also interacts with the 
surrounding environment by entrainmem, 
The parameterized entrainment rate, yp. 

is taken from experimental studies of 
buoyant jets and bubbles and is given by 

the expression 


where M is the mass of air, 2 is height, 
and R, is the radius of the core of the 
updraft. Further details of the dynamic 
and microphysical model are given in the 
above cited reference. 

{n addition to the rain cloud properties 
already mentioned, calculations of the 
concentration per unit mass of debris at 
each level in the cumulus cioud are in- 
The only 


process that alters the concentration of 


cluded in the numerical model. 


bene epee cmme 


Bet MS te 


Je om 


Se, 


den be 


debris within the sample parcel is entrain- 
ment of environmental air with a differ- 
ent debris Concentration; therefore the 
concentration per unit mass of debris in 
the cloud at level 2, Xg, Can be related 

to the concentration per unit mass of 
debris in the cloud at the previous lower 


level, Xp-1? by the expression, 


Xp. tu, dz 
Xp “Tp ae? 


where C, is the concentration per unit 
mass of debris in the environment outside 
the cloud at level %, 4 is the mass entrain- 
ment rate, and 4z is the thickness of the 
layer. At the lowest level of the model 
(at base of the cumulus cloud) the debris 
concentration is assumed to be equal to 
the environmental debris concentration, 
The numerical value of the entrainment 
rate, uw, depends on the updraft radius 

of the cloud which is, in turn, dependent 
upon the other physical properties of the 
rain cloud and the environmental sounding. 
Therefore, each cumulus cloud has a 
unique vertical distribution of debris just 
as each rain cloud has its own vertical 
distribution of drop concentrations, wp- 
draft radius, vertical velocity, etc. The 
question of whether debris is collected by 
drops or droplets is not considered in 
this simple model although it is of crucial 
importance when considering the scaveng- 
ing of debris. 


DESCRIPTION OF THE PRECIPITATING 
CLOUD 

The radiosonde sounding for Salt Lake 
City on 16 May 1971 at 1200Z was used to 
describe the environment in which pre- 


cipitating clouds developed. Scattered 
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showers were occurring in the area at 
this time. The radar-measured tops of 
rain clouds were between 6000 m and 
9000 m, and their bases were at 2700 m. 
The assumed updraft radius at cloud base 
was ! km and the assumed updraft velocity 
at cloud base was 2 m sec, The 
freezing temperature for all drops and 
droplets was assumed to be -25°C, a 
temperature which occurred at an 
altitude of about 6100 m. The thick- 
ness of the layers in the model was 
200 m. F 

The environmental and in-cloud tem- 
perature profiles are displayed in Fig. 3. 
At low levels the in-cloud temperature 
is higher than the environmental tempera- 
ture, primarily because of the release 
of latent heat as water vapor condenses, 
This positive buoyancy accelerates the 
parcel upwards as can be seen in the plot 
of vertical velocity in Fig. 4. The jog 
in the in-cloud temperature profile at 
6100 m is due to the sudden freezing of 
all drops and droplets at this level and 
the resultant release of the latent heat of 
fusion. This positive buoyancy leads to 
a corresponding increase in vertical 
velocity ai levels just above 6109 m and 
indicates the importance of freezing in 
the dynamics of cumulus convection. 
Above 7300 m the buoyancy of the parcel 
is negative and the vertical velocity drops 
to zero at 7900 m, which is, by definition, 
the top of the rain cloud. Two other 
factors also tend to decrease vertical 
velocity: the drag of drops and droplets 
which are carried along with the parcel, 
and entrainment of stationary environ- 
mental air. A cloud top of 7900 m is 
consistent with radar measurements of 
the heights of rain clouds in the Salt Lake 


———— Cloud top 


Cloud 


Height — km 


rN 


3 
Cloud base ~————_m» 
2 i { L ee (cere 
230 240 250 260 270 280 
Temperature — °K 
Fig. 3. Variation of temperature with height in the environment and in the model 


cloud, 


OU pe 8 gS yn ot 1 
2 4 6 8 10 


Vertical velocity -—— m/sec 


Fig. 4. Distribution of vertical velocity with height in the model cloud, 
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City area ai the time of the radiosonde 
sounding. 

The vertical distributions of cloud 

droplet mixing ratio and raindrop mixing 
ratio of hydrometeor water are shown in 
Fig. 5. Cloud droplets are formed by con- 
densation of water vapor in the rain cloud 
and are removed by collection by hydro- 
meteors. Raindrops are formed from 
cloud droplets and grow by coalescence. 
There is no mechanism in the model for 
a decrease of hydrometeor water in the 
sample parcel; therefore the raindrop 
water content must increase with height 
and be a maximum at the top of the cloud. 
In nataral rain, fallout of drops would 
greatly alter this particular profile. The 
discontinuity in the slope of the cloud 
droplet concentration curve at 6100 m is 
due to the freezing of drops at that level, 
and to changes in certain microphys- 
ical parameters at the freezing level. 
If all the hydrometeor water in the 
rain cloud at any instant fell to the 
ground, there would be 1.6 mm of 
rainfall. 

The characteristics of cumulus clouds 
that affect the vertical distribution of 
debris most directly are the updraft 
radius and the entrainment rate; these 
The 
mixing rate is greatest in the region of 
smallest updraft radius which corresponds 
to the height of the largest updrafts, At 
levels where the vertical velocity is 
small the entrainment rate is also small. 


quantities are plotted in Fig. 6. 


At the top of the rain cloud approx- 
imately two-thirds of the mass in the 
Sample parcel is air that has been 
entrained by the rising parcel; the 
other one-third is air that entered at 


cloud base. 
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DEBRIS DISTRiBUTIONS 


The environmental debris distribution 
that was used for the model calculations 
was taken from results of a 2BPUFF 
diffusion code (Crawford, 1966)° ealcula- 
tion. The debris concentration 100 km 
downwind from the detonation in the free 
atmosphere of a 1~-kt all-fission device 
has a Gaussian profile with a peak con- 
centration of 5107? Ci/m, The height 
of this peak concentration is 2800 m and 
the standard deviation of the profile in 
the vertical is 500 m. It is assumed that 
the variations of concentration in the 
horizontal can be neglected over the 
region that interacts with the model 
cumulus cloud, and values at the center 
of the debris distribution in the horizontal 
are used, 

The vertical profiles of debris con- 
centration in the environment and in the 
rain cloud for a 1-kt free-air burst which 
interacts with the model cloud after 
moving downwind 100 km are shown in 
Fig. 7. 


debris concentration occurs at the base 


Since .he peak environmental 


of the rain cloud, the debris concentration 
in the cloud is large. Mixing of cleaner 
environmental air as the parcel rises 

¢ reases the debris concentration in the 
cloud, but because entrainment is a 
relatively inefficient process the debris 
concentration remains high at all levels. 
Since the cumulus cloud model does 
nothing to alter the debris concentration 
below the rain cloud base, both profiles 
are identical below 2700 m. The enhance- 
ment factor, which has been defined as 
the ratio of the vertical integrals of the 
two profiles in Fig. 7, has a value of 
3,04, This implies that if all the debris 


Height — km 


Q 0.4 9.8 1.2 1.6 2.0 
Mixing ratio — g/kg 
Fig. 5. Vertical distribution of liquid water components (cloud droplets and rain drops) 
in the model cloud. 


Entrainment factor 
90.12 0.16 0.20 0.24 0.28 9.32 


Updroft radius 
(scale below) 


£6 
I 
wr Entroinment factor 
2 5 (scale above) 
x 

4 

3 

2 


9 0.4 0.8 1,2 \.6 2.0 
Updraft rodius — km 


Fig. 6. Variation of entrainment factor and size of the updraft with height for the 
model cloud. 
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Neight — km 


Environment 


0 1 2 


3 4 5 


Debris concentration — (x1073 Ci/m*) 


Fig. 7. Vertical distribution of debris in the environment and in the model cloud for 


a 1-kt all-fission free-air burst, 


The debris is assumed to interact with the 


rain cloud after moving downwind 100 km with the rain cloud where the height 
of the peak concentration for the Gaussian environmental profile is at 2.8 km 
and the standard deviation in the vertical is 0,5 km, 


in and below this cumulus cloud were 
brought scavenged, approximately three 
times the vertical integral of debris in 
the environment outside the c, mulus 
cloud would be brought to the ground with 
the rain. 

It is apparent from the previous example 
that the d oris concentration at the base 
of the r‘ in cloud is a very iniportant factor 
in controlling the debris concentration in 
the entire rain cloud. Therefore, addi- 
tional sample calculations hey been 
performeg for several other environ- 
mental debris distributions interacting 
with the same cumulus cloud; the only 
change is in the height of the peak debris 


concentration. These results are 
illustrated in Figs. 8 through 12 and are 
summarized by the data in Table 3. 

From the data it can be seen that the 
enhancement factor is largest when the 
maximum value of concentration occurs 
at the base of the sain cloud. When the 
peak concentration is below cloud base 
the enhancement factor is still greater 
than Ll, except for the limiting case when 
all the debris is below cloud base, then 
the two profiles are identical and the 
enhancement factor is 1.6. 

When the height of che peak environ- 
mental concentration is found at greater 


levels avove the base of the rain cloud, 
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Height — km 


Environment 


0 1 2 =ass ad 4 5 
Debris concentration — rodioactivity volume ; 


Fig. 8. Vertical distribution of debris in the environment and in the model cloud for 
a hypothetical free-air burst with the assumed height of the peak concentra- 
tion of the Gaussian écnvironmental profile at 2.2 km and the standard devia- 
tion in the vertical of 0.5 kin. Numerical values for the ordinate would depend 
on the total amount of radioactivity released by the free-air burst. 
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Environment 


0 i} 2 3 4 5 
Debris concentratian — radioactivity /volume 


Fig. 9. Same aa Fig. 8, except that the assumed height of the peak environmei.tal 
concentration i8 1.6 km. 
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Height — km 


Fig. 10, 


Height — km 


Fig, 11, 


Environment 


1 2 3 4 5 


Debris concentration — radioactivity /volume 


Same as Fig. 8, except that the assumed height of the peak environmental 
concentration is 3.4 km. 


Environment 


1 2 3 4 


Debris concentration — radioactivity /volume 


Same as Fig, 8, except that the assumed height of the peak environmental 
concentration is 4.4 km. 
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Debris concentration —~ radioactivity /volume 


Fig. 12. Same as Fig. 8, except that the assumed height of the peak environmental 


concentration is 5.0 km. 


Table 3, Important parameters for comparison of debris profiles, 


Standard deviation 


Height of peak in the vertical of Concentration 
environmental environmental at cloud base 
: concentration concentration 3... +3 Enhancement 
Figure (m) (m) X10" Cim “) factor 
2800 500 4,90 3.04 
8 2200 500 3.03 2,32 
9 1600 500 0.44 1.20 
10 3400 500 1.88 1.60 
11 4400 500 0.02 0.73 
12 5000 500 0 0.60 
13 4100 1000 1,88 1,21 
14 4500 1000 0,99 0.88 
transport of debris into the cumulus cell enters through cloud base, At levels 


by entrainment becomes more important, 


since only a small amount of debris 
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where the debris concentration in the 


environment is greater than the debris 


concentraticn in the rain cloud, entrain- 
ment increases the debris concentration 
in the cloud. Conversely, at levels where 
the debris concentration in the environ- 
ment is smaller than the debris concen- 
tration in the rain cloud, entrainment 
decreases the concentration of debris in 
the cloud. Consequently, the peak con- 
centration in the precipitating cell occurs 
above its base and also above the height 
of the peak environmental concentration. 
Also, the concentration of debris in the 
rain cloud in lower regions is smaller 
than the concentration of debris in the 
environment, but at higher levels, once 
again, the precipitation region contains 
more debris, The value of the enhance- 


ment factor, which is maximum when the 


9 


8 Environment 


Height — km 


0 1 2 
Debris concentration — radioactivity /volume 


Fig. 13. 
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environmental debris concentration has 
its peak value at the base of the rain cloud, 
decreases as the height of the peak 
environmental concentration increases, 
The minimum value of the enhancement 
factor of 0 when the debris is completely 
above the raincloud. The enhancement 
factor is greater than | only when signif- 
icant amounts of debris enter the pre- 
cipitating system through its base. 
Additional representative calculations 
of debris profiles for the same model 
cumulus cloud are illustrated in Figs. 13 
and 14; here the standard deviation of the 
environmental profile has been doubled, 
These results can be compared with two 
cases of the previous results, one for 


which the concentrations at the base of the 


3 4 5 


Same as Fig. 8, except that the assumed height of the peak environmental 
profile is 4.5 km and the Gaussian standard deviation is 1.0 km. 


Environment 


Height — km 


0 I 2 


3 4 5 


Debris concentration — radioactivity /volume 


Fig. 14. Same as Fig. 8, except that the assumed height of the peak environmental 
profile is 4,1 km and the Gaussian standard deviation is 1.0 km, 


rain cloud are the same (Figs. 10 and 13) 
in which case the enhancement factor is 
smaller for the broader profile, and the 
other for which the peak concentrations 
are at the same height (Figs. 11 and 14), 
in which case the enhancement factor is 
larger for the broader profile. This 
latter result follows primarily from the 
fact that the value of concentration at the 
base of the precipitation system is much 
larger for the broader environmental 
profile. 

In summary, it has been shown that 
values for the enhancement factor will 
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typically range from 0 to values of 3 or 


4. The highest values occur when the base 


of a rain cloud forms near the level of 
peak concentration, in which case the 
amount of debris available for scavenging 
can be several times <s large as the 
amount of debris that would have been 
present if the dynamics of conve-tion had 
been ignored. When the debris is mostly 
above the base of the rain clouds the 
enhancement factor is less than | imply- 
ing that the dynamics of convection tead 
to reduce the amount of debris that can 


be scavenged. 


A Microphysical Description of Precipitation Scavenging 


Aerosol particles are removed from 
the atmosphere by different processes 
depending on their size. Particles 
larger than about 1 wm can be scavenged 
by direct impaction with raindrops, This 
inertial capture process is the primary 
removal mechanism for washout or below- 
cloud scavenging, the washout efficiency 
being a function of particle size. 

For rainecut or in-cloud scavenging, 
other processes must be considered in 
addition to inertial capture, Aerosol 
particles with perfectly wettable surfaces 
and radii greater than 0.1 um can serve 
as nucleation sites for the formation of 
cloud droplets. Soluble particles with 
radii as small as 0.01 um can be 
nucleated in clouds, According to argu- 
ments by Fletcher® and experimental re- 
sults by Twomey", an insoluble particle 
whose surface departs from being com- 
pletely wettable, i.e., the contact angle 
is greater than say 6°, canno: serve as a 
condensation nucleus for typical cloud 
supersaturations, Nucleated particles 
can be deposited on the ground either 
through growth by diffusion to raindrop 
size and falling (a very rare event) or 
growing to several microns in radius and 
being accreted by a raindrop. 

Particles of radius less than about 
0,02 um display considerable Brownian 
motion under atmospheric conditions, 
allowing them to collide and attach to drop- 
lets or raindrops. This scavenging 
mechanistn is directly effective only for 
short times after the particles are formed 
since, given time, they will readily attach 


to larger debris particles or to natural 
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aerosol particles, Those partcles 
attached to droplets due to Brownian 
capture can be removed to the ground 
through accretion of the droplet. 

Fallout particles from free air bursts 
consist mainly of oxides of the dorninant 
casing materials, The various radio- 
nuclides make up a small percent of the 
totai particle mass. Indications are th it 
particles are insoluble in water and are 
wettable due to their large surface energy 
comparedtothat of water. Some of the 
radionuclides are sc’uble although it is 
not clear if the amount of soluble material 
on the particle surface can significantly 
change the nucleation characteristics as 
suggested by Hicks.?2 

Indications are that the particle size 
resulting from free air bursts are distrib- 
uted log-normally and that the mean particle 
size decreases with increasing yiela,2%44 
It appears reasonable to expect low yield 
air bursts to produce particles suscep- 
tible to inertial capture and nucleation 
scavenging. However, higher yields may 
produce particles below the nucleation 
threshold, leaving Brownian scavenging as 
the dominant removal mechanism, 

In the remainder of this section an 
attempt is made to formulate th? dominant 
scavenging processes discussed above, 
The model will be concerned with 
describing the removal of debris from a 
parcel of air in which the rain rate is 
constant throughout, it is further assumed 
that the concentration of debris particles, 
cloud droplets, and raindrops can be con- 
sidered spati: lly uniform inside this 


parcel, 


Each of the dominant removal mechanisms 
involve the collision and subsequent col- 
lection of a debris particle, or a water 
droplet which contains the particle, by a 
raindrop, Denoting as n(r, t) dr the 
number of these target spheres of radius 
between r and r+dr per unit volume of air 
at time t, the rate at which such spheres 


are removed by falling raindrops is: 
an(r, t) _ 
ome) = <a dr, t) nr, 2) q) 


A(r, t) is the tate at which one target 
particle of radius r 1s removed from the 
volume at time t, A solution to this 
equation can formally be written as: 


t 
n(r,t) = n(r, 0) exp {-f dt! Ar, of, (2) 
(eo) 


where 


n(r, 0) is n(r, t) evaluated at t = 0, 


An expression for A(r, t) as given by 


Chamberlain! can be written as: 


fa) 


Atr, t) f dR N(R, t) 7R2(V(R) - virll 
° 
X E(R, r,t). (3) 


Here N(R, t) is the number of raindrops 
per unit volume with radius R and fall 
velocity V(R). E(R, r, &) is the collision 
efficiency of the raindrop with a target 
sphere of radius r and mass density ¢&. 
vir) is the fall velocity of the target sphere, 
Scavenging by diffusiophoresis and 
thermophoresis does not appear to be 
importaut, The thermophoretic force on 
a submicronic particle juxtaposed with a 
droplet u':dergoing diffusional growth 
dominates the diffusiophoretic force g 
causing the particle to be repelled from 
the droplet, Although the attractive dif- 
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fusionhoretic force dominates for particle 
sizeS greater than about 1 wm, the re- 
sulting increase in the collision efficiency 
when the diffusiophoretic force is 
accounted for is of little significance, +” 

The removal rate given by Eq. (3) de- 
pends on tiine through the time dependence 
of the drop spectra. This expresses the 
fact that throughout the history of a cloud 
the number of drops of radius R in any 
unit volume changes with time, Under the 
assumption that the drop size distribution 
is time independent, Eq. (2) reduces to: 


n(r, t) = n{r, 0) exp [-A(r)t]. (4) 


For computation of the removal rate, A(r), 
the terminal velocities of the drop given 
by LeClair, et ai,)8 are used for Reynolds 
number less than 400 and for larger drops 
those of Gunn and Kinzer.*® Terminal 
velocities of the collected sphere are 
assumed to obey Stoke's law and the col- 
lision efficiencies used are those of 
Langmuir,2° 


from results by Best 


The drop spectra is taken 


a1 and can be written 


Here n = 2.25, W is the mixing ratio and 
ais a constant dependent upon rain rate, 
Vaiues of W and a qucted by Best”? for 
various rain rates, I, are given in 
Table 4. The drop spectrum is defined 
here Such that N(R) dR represents the 
number of drops per unit volume with 
radius between R and R+ dR. 

Debris particles of radius greater 
than about 1 um are subject tc inertial 


Table 4. Values of a and W correspond- 
ing to different rain rates, L 
From Ref, 21. 


mm 
Lie as a(cm) WX 10° (em/em)® 

0.5 0.1107 37 
1.0 0.13 67 
2.5 0.1608 145 
5.0 0.1888 261 

10,90 0.2218 470 

25.0 0.2743 1020 


ne 


capture by raindrops, Mathematically this 
situation can be described by Eqs, (1)-(4) 
where the collected particle is taken as a 
spherical debris particle. From Eq. (4) 
the fraction of debris particles of radius r 
which have been removed from the parcei 
at time t is: 


n(r, 0) - n(r, t) 


F,(r, t) : n{r, 0. 


=1-exp(-Aj(rit]. (6) 


Assuming that every particle captured by 
a drop gets deposited on the ground, 
Fy(r, t) gives the fraction of particles of 
radius rf Originally im the air that reach 
the ground as a function of time. 

Considering inertial capture to be the 
only scavenging mechanism for washout, 
Eq. (6) gives the fraction of particles of 
radius r removed from a below-cloud 
parcel a5 a function of time, Figure 15 
gives calculated values for the washout 
efficient, A,(r), as a function of particle 
size for various rain rates assuming a 
particle mass density of 5 g/cm, These 
values substantually agree with those found 
by Chamberlain, '® 

Values of the washout coefficient are 
sensitive ta the size distribution of drops, 


Using values of the drop spectra found by 
Mason and Andrews@” for warm frontal 
showers, the washout coefficients are about 
a factor of 2 smaller than those shown in 
Fig. 15. This indicates that the washout 
coefficient depends on type of rain shower 
as well as the rain rate. 

In the present model, nucleation 
scavenging is calculated by assuming every 
debris particle of radius 0.1 m< r<10um 
serves as a condensation nucleus. The 
resulting droplet is assumed to immediately 
grow to 10 um radius, Since the mass 
density of debris particles is typically 
greater than the density of water and of 
most natural aerosol particles, and 
further since there may be a Significant 
number of debris particles with a radius 
of several microns, the density of the 
debris nucleated droplet is taken as: 


rofeystesy (22-8) f yg. 


Here & is the density of the debris parti- 
cle of radius r and Cy is tre density of 
water, The radius of the nucleated drop- 
let is rye 

The fraction of debris particles of 
radius r removed from the parcel by 
nucleation and subsequent impaction with 
a raindrop is given in analogy with Eq, (5) 
and under the above assumptions as: 


F,(r, t) = 1 - exp (-Ap{rt). (8) 


Ay(r) is the removal rate for particles of 
radius 10 um and density aS shown in 

Eq. (7), Figure 16 represents the rainout 
coefficient, Ag(r), for 0,1 um < r< 10 pm, 
For particles of radius greater than 10 ym, 
the rainout and washout coefficients are 
equal size since € = q =3 gfem®, For 
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Washout coefficient ~ sec! 
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Fig. 15. Washout coefficient versus particle radius for various rain rates, 


a given rain rate, Fig. 16 shows the rain- 


out coefficient to be relatively constant 
for the particle sizes considered. 

In order to make a complete rainout 
assessment, one should have details of 
the size distribution of debris particles. 
However, as a result of the assumption 
that every debris particle between 0,1 
and 10 um radius serves as a condensa- 
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tion nucleus when subjected to a cloud 
environment, a further simplification 
concerning the rainout calculation becomes 
apparent. Let m(r) dr represent the 

mass contained ir debris particles of 
radius between r andr+ dr. The total 
mass of debris deposited on the ground by 
particles of radius greater than 0.1 um 
becomes: 


-1 


w& 


no 


Rainovt coefficient — sec 


Rain rate = 25 mm/hr 


1.0 1. 1.4 1. 


Particle radius — (x 1073 em) 


0.5 mm/hr 
a. 
0 co ee Eeeeeen Eee 
i) 0.2 0.4 0.6 0.8 
Fig. 16. 
(3) 
M(t) = dr m(r) 
0.i um 


x [1 - exp (Ay(r)t)) ° (9) 
Here Ay is the rainout coefficient for a 
particular rain rate as plotted in Fig. 16. 
Since A,(r) is a weak function of r, Eq. (9) 
can be approximated by: 
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Rainout coefficient versus particle radius for various rain rates, 


t 00 
“AS 
Mit) = (: -e Nf dr m(r) 
0,1 
t 
-A 
a eS 2 
= ( e ) M.- 


Here M, is the mass contained in particles 
of radius greater than 0.1 4. Under the 
further assumption that the radioactivity 
is contained in particles of greater than 


(10) 


0.1 um radius, which is approximately An evaluation of the potential hazard 
true for yields as low as 1 kt, the fraction due to rainout from a 1-kt free-air burst 
of radioactivity deposited on the ground as __— baSed on the above results can now be 

a function of time is plotted in Fig, 17. made. Consider the debris system to 
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Fig. 17, Fraction of radioactivity originally in the rain cloud deposited on the ground 
versus time for various rain rates, 
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Initial infinite whole-body exposure rate — 2/hr 


ty) 


Fig. 18. 


© 25 mm/hr rain rote 


0 2.5 mm/hr rain rate 


100 1000 
Distance from ground zero —~ km 


Data points represent initial gamma exposure rates due to rainout of debris 
from a 1-kt free-air burst versus distance from ground zero, Rain is 
assumed to fall for 15 min with rain rates of 2.5 and 25 mm/hr., The solid 
curves are reproduced from Fig, 1. 
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move downwind and to diffuse according 

to the Lagrangian diffusion code (2BPUFF) 
under non-shezr conditions with horizontal 
dissipation, € = 1.8 ergs/g-sec. At vari- 
ous distances downwind the debris system 
is assumed to interact with a cumulus 
cloud, as described previously, giving an 
enhancement factor of 3, The radioactive 
particles are assumed uniformally 
distributed within the rain cloud and 
deposited by rainout as a function of time 
and for various rain rates as given in 
Fig. 17, 
assumed that the drops fall on the area 
directly below the rain cloud. The data 
points in Fig, 18 show values of the initial 


Washout is neglected and it is 


infinite whole body exposure rate versus 
distance downwind for rain rates of 25 and 
2.5 mm/hr lasting for 15 min. The curves 
in Fig, 18 are reproduced from Fig. 1. 
These new results on the potential 
hazards from rainout of debris from the 
1-kt cloud are worthy of special comment, 
namely, that the results of the present 
analysis indicate that the potential dose as 
a function of range for rainout (Fig. 18) 
constitutes a reasonable estimate of the 
expected hazard for the conditions cited. 
This expected estimate level of hazard 
is quite comparable to that previously 
reported in "Potential Exposures from 


Low- Yield Free Air Bursts," During 
the remainder of this ccntract year, the 
LLL Rainout Team will continue to refine 
the present estimate, as represented in 
Fig. 18, for effects of wind shear (not 
expected to be large in summer), and 
other assumptions inherent in the analy- 
sis. Further work is required in regard 
to the variation of particle size distribu- 
tion with yield in order to extend the 
present analysis techniques to yields other 
than 1 kt. This work will be a portion of 
our proposal for the coming year. In the 
future we will also evaluate the probability 
associated with meteorological conditions 
consistent with Fig. 18 and other scenarios 
of interest to DNA. 

Rainout estimates assuming nucleation 
scavenging and inertial capture were 
possible for a 1-kt free-air burst be- 
cause in this case it appears that very few 
of the debris particles are smaller than 
0.1 um radius. 22 It appears that the mean 
particle size decreases with increased 
yield such that for higher yields a 
Significant portion of the radioactivity is 
contained in particles smaller than0.1 ym. 
This being the case, different scavenging 
mechanisms (probably Brownian capture) 
become important and the above analysis 
does not apply. 


Thoughts on Altering the Scavenging of Debris Particles 


The problem here is to change or te 
control the rate of rainout from clouds. 
We are studying the possibilities of 
changing their macrophysical dynamics 
by interfering with their microphysical 
processes, One can add toa cloud or 
take parts of it away, and at first look 
it appears that adding something is the 


easier Operation. 
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There are a number of possible effects 
of additives on the microphysical proc- 
esses of water and ice particle formation 
in a cloud that can influence the overall 
macrophysical dynamics of the cloud. 
Additives in the form of an excess num- 
ber of nucleating agents could hinder the 
development of water and/or ice particles 
to a size large enough to precipitate out 


of a cloud, Small amounts, a monolayer 
or even a partial monolayer, of materials 
with low surface energies adsorbed on to 
the high energy surfaces of nucleating 
particles can change their wetting charac- 
teristics and thus their nucleating rates 
and possibly put them into the nonnucleat- 
ing category as far as the cloud is con- 
cerned, Small amounts of materiais can 
also affect the electrical properties of the 
water and/or ice particles, the elastic 
and inelastic collision processes, the 
accommodation and reflection coefficients 
for the water molecule on water and ice, 
and change the condensation and/or evap- 
oration rates of water molecules between 
water droplets and between water and ice 
particles, 

Let us look at the use of additives to 
modify the rate of nucleation of particles 
by changing their wetting characteristics. 
We shall assume the validity of the 
Fletcher analysis® of heterogeneous nucle- 
ation of water vapor by insoluble particles 
as a function of partinle size, contact 
angle and saturation ratio in the cloud. 
Figure 19 is a set of Fletcher's curves 
for a given nucleation rate. Since the 
saturation ratio in most clouds is about 
1,001 and only under unusual circum- 
stances may reach as high as 1,01, it can 
be seen from Fig. 19 that a particle with 
radius less than 10° cm will not act as a 
nucleation center. Even particles with a 
radius greater than 10° cm need to be 
almost completely wettable before they 
can act as nucleation centers for water 
vapor in a typical cloud, 

The solid particles that are expected 
to be picked up or injected into a cloud 
would be metal oxides in the form of crys- 
tals or glasses and possibly some free 
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metallic particles, All such particles 
have a large surface energy compared to 
See Tables 5 and 6. 


One would expect water to wet and spread 
29-31 


water (73 ergs/cm“). 
on all such particles, and from the 
Fletcher analysis all such particles would 
be active nucleating agents in a cloud if 
they were 10° cm in radius or larger.? 
The rates of nucleation by these solid 
high surface energy particles can be 
changed by the addition of low surface 
energy materials that will adsorb onto 
the surface of these particles in the cloud 
and lower their surface energies to less 


than that of water. A monolayer of 


Table 5. Surface energies of various 
materials, taken from refer- 
ences cited. 

y at 0°K 
Materials (erg/cm?) Reference 


Oxides of General Formula MO 


MgO 1090 24 
FeO 1060 24 
MnO 1010 24 
CaO 820 24 
SrO 700 24 
BaO 604 24 
BeO >1420 24 
cdo 530 24 
LnO 600 24 
PbO 250 24 


General Forraula M,03 


B,0, 79.5 @ 900°C 25 
Al,O, 905 @ 1850°C 26 
General Formula MO, 

UO, 642 + 20% 24 
ZrO, 800 + 20% 24 
ThO, 530 t 20% 24 
SiO, >800 24 
TiO, >8006 24 


Saturation ratio 


6 


1978 1077 10 


1075 10 


Radius ~ cm 


Fig. 19. Critical supersaturation ratio for nucleation of a water droplet in 1 sec on a 
particle of given radius and surface properties defined by m = cos ¢ whered¢ 


is the contact angle (from Ref. 3). 


perfluoro fatty acid on platinum reduces 
its specific surface free energy from 
about 2000 ergs/em” to about 10 ergs/ 
em?,32,33 The water should not wet nor 
spread on these particles. Even a smali 
departure from perfect wettability will 
hinder the nucleation dynamics, There 
are a number of possible materiais that 
could be selected, and these could be 
introduced into the cloud by rockets, bal- 
loons, smudge fires, antiaircraft guns, 
the device, etc. If one could make the 
additive available to be adsorbed onto the 
particles not much material would be 
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needed. From Table 7 it can be seen that 
for each em? of metal or metal oxide, 
subdivided into 10'° particles with sides 
of 10° cm, it would take only 10° mole 
of adsorbable material to give all parti- 
cles a monomolecular coat, 

For optimum control of rainout one 
may want to use one or more types of 
additives at the same time, It may be 
best to keep the insoluble particles from 
wetting by an adsorbable additive and to 
overnucleate the cloud with NaCl parti- 
cles with masses equal to or greater than 
102 grams per particle. This would 


ak) Seti 


cd 


keep the insoluble particles from nucleat- 
ing and furnish enough water soluble par- 


ticles for nucleating centers to keep the 
droplets of water and particles of ice 
small enough so that they will not precipi- 
tate out and carry out the solid particles 
by impact, 

There are a number of possible effects 
of additives on cloud dynamics about which 
little or nothing is known, Research into 
these effects is desperately needed, 

It should also be kept in mind that a 
cloud is a dynamic system, and the appli- 
cability of equilibrium arguments is some- 
what limited and may at times be mislead- 
ing. The Knudsen equation, for instance, 
gives too large a gas transport to and 
from a water drop under nonequilibrium 
conditions, It is probable that the con- 
densation and evaporation coefficients of 
water molecules at liquid water surfaces 
can be modified by as much as a factor 
of 10 by contaminating the water sur- 
faces?4~36 or even aging them. There 
are many different combinations of fac- 
tors involved in the interaction between 
incident water molecules and surfaces of 
water drops and ice particles, These 
factors are the: 


Surface area of a 1-cm cube with subdivision, 


Table 6, Surface energy of some metal- 
lic elements (at their melting 
points),¢!; 

Melting Surface 
point energy, ¥ 
Metal (°C) (ergs /cm2) 

Aluminum 660 900 

Antimony 630 370 

Beryllium 1400 1000 

Cadmium 321 620 

Cobalt 1495 1530 

Copper 1083 1100 

Gold 1063 1120 

Iron 1534 1500 

Lead 325 450 

Magnesium 650 560 

Mercury -39 460 

Nickel 1453 1700 

Platinum 1769 1800 

Silver 961 920 

Tin 232 570 

Tungsten 3410 2300 

Uranium 

Zinc 420 790 

Table 7, 
Length of side Number of 
(cm) cubes 

10° 10° 

10} 10° 

1072 10° 

1973 10° 

io? (= 1pm) 107? 

vw? 1015 

10° i918 

10°" 101 

10°8 (31 A) 


aw ,gssume that it takes 10 


6 xX 10-4 molecules, 
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Total surface 


Moles of material 


area (cm*“) to cover surface 


6 x 10° 10 19 
6x10! 10? 
6 x 102 108 
6 x 10° 107 
6 x 104 10 6 
6x 10° 10° 
6 x 108 104 
6x 10" 10°3 


2 
molecules to cover each cm” and each mole has 


specific condensation coefficient, 
specific evaporation coefficient, 
thermal accommodation coefficient, 


momentum accommodation coef- 
ficient, 

@ extent of compliance with the cosine 

law for diffuse reflection, and 

e@ extent the cosine law would obtain 

due to scattering from surface 
irregularities alone, 

Even though Young?" gave a correct 
form of the interaction of a liquid at equi- 
librium on a solid surface in 1805 it was 
only in 1959 that Johnson®® gave a clear 
statement of the problem and a thermo- 
dynamic justification of Young's equation 
along the lines of Gibbs' derivation. The 
notion of a wetting tension was given by 
Dupré (1869),°9 and Cooper and Nuttall 


(1915)?9 gave the conditions for spreading 


and nonspreading. Harkins (1920-1941)*! 
further developed the notion ~f a spread- 
ing tension more from the point of view of 
a gradient of the free energy of the sys- 
All of the 
above developments are for equilibrium 


tem with respect to the area. 


conditions and assume no electrification 
of the system, The molecular and/or 
atomic interaction mechanisms of wetting 
and wetting rates are not well understood, 
If one considers small drops on surfaces, 
ignoring inertial effects, one might expect 
the wetting rate of a surface to be ex~ 
pressed by equating the rate of change of 
the reversible work due to "capillary" 
forces to the rate of energy dissipation 
due to the viscous flow of the film over 
the surface. There are a number of efforts 
needed to understand the off-equilibrium 
dynamics in the microphysics of clouds. 


Dry Deposition 


One method by which radioactivity 
from air burst debris clouds may be de- 
posited upon the ground, generally over 


long periods of time, is through the process 


of dry deposition, In this process the very 
small particles inthe clouds are dispersed 
downward through the atmosphere by dif- 
fusive eddies. As the small particles 
reach the ground, they become attached to 
various surfaces which they may encoun- 
ter and, unless removed, become sources 
of radicactivity at these points. The 
amount deposited upon a given surface 
depends upon a so-called velocity of dep- 
osition (or deposition velocity) defined? 5 
ag 


Amount deposited per em? 


Vo = of surface per sec 
c olumetric concentration 


er cm’ above surface 
p 
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where the volumetric concentration is 
measured very near the surface. Values 
of V, vary from 0,1 to 10 em/sec with 

most reported values near 1 cm/sec (see 
Ref. 15), 


deposition, D, one needs a time history of 


Thus, to determine the total 


the volumetric corcentration x near the 
surface in order to evaluate the integral 


ty 
x V,dt. 


Bal 
Jj 


is 


Of course decay must be accounted for 
when considering radioactive materials. 

When the debris cloud from a free-air 
burst becomes stabilized and moves down- 
stream, the concentration of radioactivity 
near the ground will usually be small 


(depending upon height of stabilization) 
until the radioactive particles are dis- 
persed downward by the action of vertical 
diffusion, Since this process can be ex- 
pected to take some time and the amounts 
deposited are small for low concentra- 
tions, the potential threat from one or 
even a few bursts is seen to be quite 
small in most cases. An exception would 
be when the surface rises rapidly toward 
the higher concentrations near cloud cen- 
ter as in a ange of mountains or high 
hills, 
low-yield bursts were detonated within a 
limited area, the resulting cloud would be 


If, however, a large number of 


large in horizontal extent but would be 
The radio- 
active particles would soon become dis~ 
persed to the surface, and, due to the 
large horizontal size of the cloud, high 
surface concentrations could be present 
at a given point for long periods depend- 
ing upon the time of pessage of the cloud, 
Also, as the cloud aged and became very 
large horizontally, even a small surface 


concentration, when integrated over a 


stabilized at a low altitude, 


long period of time may result in large 
values of dry deposition; this effect would 
be further enhanced at points far down- 


stream if conditions were such that no or 
very little deposition of radioactivity 
occurred earlier, 

The effect of wind velocity, i.e., the 
mean speed at which the debris cloud is 
carried downstream, will also have a 
bearing upon the total deposition since, 
for a given distance from ground zero, 
the mean speed will determine the amount 
of decay, and to some extent, the amount 
of dispersion of the radioactivity within 
the cloud, 

Table 8 presents some calculations 
which demonstrate the effects of stabil- 
izaticr height and wind speed upon the 
infinite dose from gross gamma due to 
dry deposition over flat ground at various 
distances downstream from ground zero, 
It is assumed that the clouds are created 
by the detonation of 10 Mt of low yield 
devices within an area bounded by a 
radius of 80 km. The calculations were 
made using the computer code 2BPUFF® 
programmed with a deposition velocity of 
l cm/sec. Trials 1 through 3 demon- 
strate the relatively small exposures 
obtained from debris clouds which stabi- 
lize at heights of 5 to 6 km, The 
higher vajues in Trial 2 are caused 


Table 8. Effects of stabilization heighi and wind speed on the infinite gross gamma 
dose due to dry deposition, 
Cloud 
center Wind Initial 
Trial height speed thickness Dry deposition, R 
No. (km) (m/sec) (km) 200 km 2000 km 20,000 km 

1 6 20 4.4 neg. neg. 0,003 
2 5 20 4.4 0,002 0.8 0.2 
3 5 10 4.4 0.01 0.005 0.01 
4 2 20 4.4 40 2 0.15 
5 1 10 4.4 100 2 0.15 
6 3 10 3.0 neg. neg. 10 
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by experimental use of aigher values 
of the vertice! diffusivity. In Trial 4, 
the cloud center is lowered to 2 km and 
the thickness” is such that a reasoi.ably 
high concentration is present at the sur- 


faceatthetime of stabilization. The circum- 


stances are sufficient to cause large expo~ 
sures close-in with decreasing but high 
exposures farther downstream. When the 
cloud center is lowered to 1 km (Trial 5), 
exposures are higher at all distances 
considered even though the windspeed is 
halved and the travel time therefore 
doupled with respect to Trial 4. Trial 6 
was made to demonstrate the effect of a 
debris cloud in which diffusion and depo- 
sition are inhibited until the cloud is 
20,000 km downstream, at which time a 
large vertical diffusivity is allewed to 
bring debris to the surface. This situa- 
tion might represent the cloud's passage 


over a large cold dome of air and/or pass- 
ing over an area in which a low-level 
temperature inversion is present for a 
long period of time. Furthermore, the 
increase at 20,000 km is partly due 
to the fact that the terrain was pro- 
grammed to siope upward toward, but 
not to reach, cloud center, thereby 
simulating the effect of a debris cloud 
moving over a range of high hills or 
low mountains. 

More work needs to be done on the 
effects of dry deposition, Deposition 
velovities over different surfaces and 
velocities of different particle sizes need 
to be considered, and the inhibitory 
effects of inversions and the enhancement 
effect of strong convection and overturn- 
ing, as well az the extent to which sloping 
land increases or decreases deposition 
should be studied. 


A Three-Dimensional Atmospheric Diffusion 
Particle-In-Cell Code (ADPIC) 


ADPIC is a numerical three-dimensional 
Cartesian particle diffusion code, capable 
of calculating the time dependent distribu- 
tion of air pollutants under many condi- 
tiens, which include strongly distorted 
advection wind fields and calm conditions 
for which Gaussian models are not suit- 
able, Basically the code solves the three- 
dimensional advection-diffusion equation 
in its conservative form (pseudo-velocity 
technique), for a given mass-consistent 
advection field, by finite difference 


"The thickness of a debris cloud in which 
the debris is nermally distributed in the 
vertical, as is the case here, is the dis- 
tance between the two standard deviation 
concentration points; i.e., thickness = 40,. 
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approximations in Cartesian coordinates, 
The method is based on the particle-in- 


42-44 vith the pollutant 


cell technique 
concentration represented statistically by 
imbedded Lagrangian marker particles in 
an Eulerian grid. The inclusion of 
sources, sinks and gravity terms is 
planned, The capability of the grid to 
automatica!ly expand and travel with a 
pollutant cloud, makes the code particu- 
larly suitable for the study of single puff 
releases, such as a free-air burst. 


THE METHOP 


The pseudo-velocity method consists 
of the following: Given the nonlinear 


transport-diffusion equation, 


3Q 


at (11) 


+U.-0Q=%9. (KYQ) 
where Q is a scalar concentration, K the 
diffusion coefficient and U the (given) 
mass consistent wind advection field, we 
can, under the assumption of incompress- 
ibility, replace the U - VQ term by 

V+ (QU). Upon combining the two diver- 
gence terms we can rewrite Eq. (11) in its 
conservative (pseudo-velocity) form, 


a = Ks _ aQ 
H+ ¥- [a - §%)]-F 

+7- (QvUy=0 (12) 
where U' = T - K/Q VQ is the pseudo 


transport velocity. 

The grid mesh of the cade is repre- 
sented by an Eulerian grid consisting of 
three-dimensional rectangular cells of 
uniform size, The concentrations Q are 
defined at the centers of the cells and the 
velocities U, U’ and -K/Q VQ are defined 
at the cell corners, The locations of the 
particles, which represent the pollutant 
cloud statistically, are defined by their 
individual Lagrangian coordinates within 
the Eulerian fixed grid. 

A time cycle of the code is divided 
into an Eulerian step and « Lagrangian 
step and proceeds as follows: 

1, Eulerian Step: The concentrations, Q, 
given for each cell at the beginning of 
the cycle, are used to calculate the 
diffusivity velocities 0, = -K/Q VQ 
which are then added to the given wind 
advection velocities U to yield a pseudo 


velocity U' for each cell corner. 


Lagrangian Step: Each marker parti- 
cle contained in a given cell is trans- 
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ported for one time step At witha 
velocity U', which is computed from 
the pseudo velocities U' at the corners 
of the cell, and the particle coordi- 
nates X by a volume weighting scheme, 

X (new) = X (old) + u, At, (13) 
Finally, a new concentration distribu- 
tion, Q, is calculated from the new 
particle positions, which ends the 
cycle, 

The chief advantage of using sucha 
hybrid Eulerian-Lagrangian scheme is 
that the fictitious diffusion inherent ina 
purely Eulerian scheme is eliminated. 
The truncation errors inherent in the 
finite difference algoritnms remain, of 
course, and must be minimized by the 
choice of the time-step. 

The bounrlary conditions are imple- 
mented in a layer of boundary cells that 
surround the grid, {nm the case of a grow- 
ing puff cloud the Eulerian grid will be 
expanded automatically, once a particle 
has moved into any of the boundary cells. 
Topography and inversion layer bounda- 
ries are represented by a deposition and 
reflection boundary conditions, respec~ 
tively. The use of relative coordinates 
allows the grid to follow a pollutant cloud 
along tts trajectory, The computer time 
required for a 15 X 15 X 15 = 2375 celi 
Eulerian grid with 4000 particles is 
approximately one second per time cycle 
on a CDC 7600, 


“Reference tc a company or product 
name does not imply approval or recom- 
mendation of the product by the University 
of California or the U.S, Atomic Energy 
Commission to the exclusion of athers 
that may be suitable. 


Future Work 


Future plans include continued re- 
search into the areas outlined in this 
report. To account for shear cffects on 
the debris system, a three-dimensional 
particle in cell diffusion code capable of 
calculating the time dependent distribu- 
tion of debris under distorted advection 
wind fields is being developed and 
validated. 

The dynamics of interaction between a 
debris system and a precipitating system 
determine the amount of debris available 
for scavenging and the time history of the 
wet deposition on the earth's surface, 
Future plans include adapting existing 
one- ani two-dimensional tinie~dependent 
rain cloud models to describe the interac- 
tion of a precipitating system with a 
debris cloud. Incorporation of a micro- 
physical description of scavenging into 
these numerical] rain cloud models wil} 
enable estimates to be made of the fate of 
debris within the cloud as a function of 
time and of the spatial and temporal dis~ 
tribution of radioactivity deposited on the 
ground, 

Theoretical investigations into the for- 
mation of debris particles to determine 
the particle size distribution for various 
yields and its evolution with time are 
planned, A continued effort wil) be made 
toward understanding the chemical compo- 
sition of the debris particles and manipu- 
lation of the wettability of the particle 
surfaces to affect nucleation, 

Additional work is needed to more 
realistically describe rainout. Among 
the problems to be studied are what frac- 
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tion of debris particles can serve as con- 
densation nuclei in a given precipitation 
system and how these particles compete 
with the natural condensation nuclei, Ifa 
large number of small particles result 
from a detonation, they would attach 
themselves to natural aerosol particles 

by virtue of their Brownian motion, and 
the removal of the contaminated natural 
aerosol particles must then be considered. 

The calculation of dry deposition can 
be improved by including the effect of 
wind shear on the traveling debris system. 
Another problem is the efficiency by 
which these particles attach to the earth's 
surface, Future work into this problem 
area is planned as well as an assessment 
of the global effects of wet and dry 
deposition, 

Fohl and Zalay?® suggeSted that the 
height of rise of two buoyant fireballs 
released simultaneously and in close 
Proximity will be decreased by their 
mutual interaction. Further and more 
fundamental work is needed to evaluate 
and predict this phenomenon and to assess 
its implications on the scavenging of 
debris, 

In order to more effectively assess the 
consequences of the rainout proolem, a 
climatological study for areas of interest 
is planned. Such information as the depth 
of convection cell size, the fraction of 
time that precipitating systems are pres- 
ent, and the type of flow above precipitat- 
ing systems coupled with improved rain- 
out estimates, can be useful for military 
planning, 
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